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Abstract. Here we apply a recently developed maximum likelihood method for determining best-fit parameterised 
lens models to observations of the rich lensing cluster Abell 1689. The observations that we use were taken with 
the ESO/MPG Wide Field Imager. The wide field-of-view enables us to use the weakly lensed images of faint 
background objects on an unsurpassed range of scales, 0.12 hr 1 Mpc < R < 1.8 h~ 1 Mpc from the cluster centre, 
to determine best-fit models for the 1-parameter singular isothermal sphere (SIS), 2-parameter general power-law 
and NFW models, and 3-parameter singular isothermal ellipsoid (SIE). The best-fit SIS has an Einstein radius 
0 E = 0'.37 (0.043 /G 1 Mpc) i.e. a velocity dispersion itid = IO 28 I 42 kms -1 in an = 1.0, A = 0.0 cosmology. For 
the best-fit NFW profile, the virial radius r 2 oo = 1.14 A -1 Mpc and the concentration parameter c = 4.7, giving 
a virial mass M 200 = 5.7 x 10 14 ft -1 Mg. At q = 0.88, the slope of the best-fit power-law model is slightly flatter 
than isothermal (q = 1 . 0 ), indicating that the galaxies most important in the fitting procedure lie inside the scale 
radius r a . By fitting an SIE, the deviation of the projected mass distribution from circular symmetry is evident, 
with a best-fit axial ratio / = 0.74. 
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1. Introduction 

The gravitational field of a cluster acting as a weak lens 
distorts the shapes of background galaxies (shear effect), 
and changes their observed number density (magnification 
effect); see Mellier (1999) and Bartelmann & Schneider 
(2001) for recent reviews. These imprints on the back¬ 
ground galaxy population can be used to constrain the 
projected mass distribution of the cluster. Parameterised 
cluster models are essential when statistical comparisons 
of clusters are to be made, and their error properties are 
well understood in contrast to those of non-parametric 
mass maps. 

In Schneider et al. (2000; hereafter SKE) we devel¬ 
oped likelihood techniques to obtain parameterised models 
from weak lensing data, and quantified the accuracy with 
which parameters can be recovered using the magnifica¬ 
tion and/or shear information. This work concentrated on 
power-law models, and was extended in King & Schneider 
(2001; hereafter KS) to encompass NFW (Navarro et al. 
1996; 1997) and ellipsoidal models (Kormann et al. 1994), 
also investigating the effect of uncertainty in the redshifts 
of potential source galaxies. In King et al. (2001; hereafter 
KSS) we considered the influence of substructure, conclud¬ 
ing that the method is robust against the presence of sub¬ 
structure at a level consistent with N-body simulations. 


Gray et al. (2000) applied the magnification maximum 
likelihood method to near-infrared CIRSI observations of 
the cluster Abell 2219, and extended it to include the ef¬ 
fect of incompleteness. 

In this paper we obtain parameterised models of a 
lensing cluster by applying the shear maximum likelihood 
method. We use WFI observations of the rich z<j = 0.18 
lensing cluster Abell 1689, the observational details of 
which are presented in Clowe & Schneider (2001). 

The families of lens models considered are the singu¬ 
lar isothermal sphere (SIS), the singular isothermal el¬ 
lipsoid (SIE), a general power-law model and the NFW 
profile. The SIS model is frequently used as the basis 
model for lensing studies, and the introduction of an el- 
lipticity, yielding the SIE model, allows us to investigate 
non-spherically symmetric dark matter distributions. The 
motivation for considering the NFW profile is that it is a 
good description of the radial density profiles of virialised 
dark matter halos formed in cosmological simulations of 
hierarchical clustering. We also investigate a general sin¬ 
gle slope power-law model, since in KSS we found that 
lensing simulations through an N-body cluster were well 
fit by this model. 

The structure of this work is as follows: in Sect. 2 we 
summarise how the shear information is incorporated into 
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a likelihood analysis, and outline the relevant observa¬ 
tional parameters. Our results are presented in Sect. 3, 
and we finish with a discussion and conclusions in Sect. 4. 

2. An outline of the shear method and 
observational details 

The notation and basic relationships used in this paper are 
as in KS, and the reader is referred to SKE for details of 
the maximum likelihood treatment. Throughout, standard 
lensing notation is used (e.g. Schneider, Ehlers & Falco 
1992; Bartelmann & Schneider 2001). 

2.1. The shear method 


then the redshift distribution of the galaxies would be¬ 
come important and this sheet approximation is no longer 
robust. We take the background sources to be at z = 1.0 
unless otherwise stated. This is motivated by Clowe & 
Schneider (2001), who used the HDF-South photometric 
redshift catalogue from Fontana et al. (1999) and found 
a mean galaxy redshift of 1.15 after applying the same 
magnitude cuts as above. Since this value is affected by 
Poisson noise and cosmic variance, being measured from 
only 48 galaxies, we consider the influence of the assumed 
source redshift in Sect. 3.4, Unless otherwise stated 


we 


give the results for an Q = 1.0, A = 0.0 cosmology, for a 
Hubble constant Ho = 100/ikms -1 Mpc -1 . At z& = 0.18, 
an angular scale of l'.O corresponds to a linear scale of 
0.12 ft -1 Mpc. 


The basis of using the distorted images of background 
galaxies to constrain the cluster model is that the expec¬ 
tation value for the lensed ellipticity (e) = g , the reduced 
shear, in the non-critical regime (e.g. Schramm & Kayser 
1995), and that (e) = 1 /q* in the critical regime (Seitz & 
Schneider 1997). 

A parameterised model specifies the dependence of 
g(fl) on position d for a profile family; these are given 
in Sect. 2.2 below, for the families considered. Minimising 
the shear log-likelihood function (see SKE for details) 


iV 7 

^7 = -^2^Pe{£i\g(^i)) ( 1 ) 

i=1 

over the Ab, galaxy images gives 7 r max , the parameters 
most consistent with the probability distribution of lensed 
ellipticities p e (e| g(d)). The trial parameters en route to 
the best-fit will be referred to by n. Noise in this method 
arises from the intrinsic dispersion in the galaxy ellipticity 
distribution, a e s, and observational effects. 


2.3. Profiles 

2.3.1. The isothermal profile 

The singular isothermal profile, or SIS, is characterised by 
a velocity dispersion <tid that is constant with radius. It 
is parameterised by an Einstein radius 

= (2) 

which marks the boundary between the strong and weak 
lensing regimes. Note that for a given ctid, the Einstein 
radius depends on Q and A through the angular diameter 
distances, but that it is independent of H 0 . The surface 
mass density is given by 

«( 0 ) = 0-5 (J-^j . ( 3 ) 

The mean surface mass density inside 9 is R(9) = 2 k(9), 
and further |y| = n. 


2.2. Observational parameters 

The catalogue of lensed galaxies in the cluster A1689 con¬ 
taining positions fli and lensed ellipticities n is described 
in Clowe & Schneider (2001). Objects therein were se¬ 
lected to have 23 < rriR < 25.5, and a signal-to-noise 
greater than 9. The data field used in the method has inner 
and outer radii of 6 ln = l'.O (0.12 /i -1 Mpc) and 0 O ut = 15'.4 
(1.8 A -1 Mpc) respectively. Within this aperture, there are 
AL = 19400 images for which an ellipticity can be mea¬ 
sured, and can therefore be used for the shear method. 
These galaxies have a Gaussian probability distribution 
with 2-D dispersion cr £ = 0.424; since \g\ ^ 0.2, correction 
to the unlensed dispersion is not necessary. 

The redshift distribution of the source galaxy popula¬ 
tion can safely be neglected, since the lens is at a fairly 
low redshift (Zd = 0.18) (see for example Bartelmann & 
Schneider 2001). This amounts to approximating the back¬ 
ground galaxies to be located at a redshift corresponding 
to the mean value of their lensing effectiveness parameter 
w(z). If the cluster was at a higher redshift, *5 0.25 say, 


2.3.2. The singular isothermal ellipsoidal profile 


The SIE profile is characterised by a velocity disper¬ 
sion that is constant with radius; it is parameterised by 
an equivalent angular Einstein radius 9 e, an axial ratio 
0 < / < 1 and an orientation a. Let the polar coordinates 
in the lens plane be = ( 9cos(f> , 6 *sin</>). The equivalent 
angular Einstein radius for a lens with velocity dispersion 
<tid is given by Eq.|j. The dimensionless surface mass den¬ 
sity is given by 


k(6»,0) 

b 


V7 

2b 

-jL- ((1 + f 2 ) + (1 - / 2 ) cos[2(^ - a )]) 1 (4) 


and the magnification is 


p{0,4>) 


1 

1 - 2k(9, </>) ‘ 


(5) 


The components of the shear are 


7 i = —k cos( 2 (f>); 72 = — k sin( 2 <)>) . 


( 6 ) 
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2.3.3. The power-law profile 


The power-law profile can be characterised by an Einstein 
radius 0e and a slope q , such that 



and 



The shear is 


7 ( 0 ) 


g 

2 



(7) 

( 8 ) 

(9) 


In KSS we found that the power-law profile can provide 
almost as good a fit as the NFW profile, depending on the 
size of the cluster’s scale radius r s relative to the available 
data field. 


2.3.4. The NFW profile 


The properties of the NFW profile in the context of gravi¬ 
tational lensing have been discussed by several authors in¬ 
cluding Bartelmann (1996) and Wright & Brainerd (2000). 
We can parameterise this profile with a virial radius r 2 oo i 
and a dimensionless concentration parameter c, which are 
related through a scale radius r s = r 2 oo /c. Inside r 2 oo, the 
mass density of the halo equals 200 p c , where p c = 
is the critical density of the Universe at the redshift of 
the halo. The characteristic overdensity of the halo, 4 C , is 
related to c through 


200 c 3 

3 ln(l + c) + c/(l + c) 


( 10 ) 


Then the density profile is 


p(r) 


4c Pc 

(r/r s )(l +r/r s ) 2 ’ 


( 11 ) 


which is shallower than isothermal (r -2 in 3-D) near the 
halo centre and steeper than isothermal for r *5 r s . For 
details of quantities required in the lensing analysis, see 
KS. 


3. Results 

3.1. The best-fit parameterised models 

The shear method was applied to the catalogue of back¬ 
ground objects, finding the best-fit model parameters, for 
each of the model families described above. 

For the 1-parameter SIS, the best-fit #e=0'.37; this cor¬ 
responds to a linear size of 0.043 h _1 Mpc at the lens red- 
shift. The values of 2A£ = 2 [£{ir) — T( 7 r max )] were calcu¬ 
lated from the data in order to obtain confidence levels for 
the fit. The upper (lower) 1-cr limits are 01395 (0'.34) and 
the upper (lower) 3 -a limits are 0144 (0129). This trans¬ 
lates into trio = 1028 j( 42 kms _1 ( 1 -cr errors), for fl = 1.0 


Table 1 . The best-fit parameters 7 r max for the SIE, 
power-law (POW) and NFW profiles are presented in 
the middle column. The final column shows 2A£ = 
m ^"max (Model)) — £ ( 

T^max (ModelNFw))]- 


Model 

Best-fit parameters 7r max 

2A£ 

SIE 

POW 

NFW 

0e = 0137; / = 0.74; a = -27° 

9e = 0'.3; q = 0.88 

^200 = 1.14 h -1 Mpc; c = 4.7 

2.5 

3.2 


and A = 0.0. In an U = 0.3, A = 0.7 Universe we obtain 
aiu = 998l 33 kms _1 (1-cr errors); this is within the 1-cr er¬ 
ror of the result given in Clowe & Schneider (2001) for the 
same choice of cosmology and background source redshift, 
calculated with a direct % 2 fit to the shear profile. 

In Table [j], we present the values of 7 r max most 
consistent with the data for the 2 -parameter power- 
law and NFW profiles, and for the 3-parameter SIE 
profile. Comparing the 2-paranreter families, the over¬ 
all best-fit model comes from the NFW family but 
it is only marginally preferred over the power-law 
model. For 2 degrees of freedom, the values of 2A£ = 
2 {£ ( 7 r max (power — law)) — £ (7r max (NFW))] show that the 
best-fit models are only distinguished at about 80% con¬ 
fidence. 

In order to fit the SIE profile, we adopted a complex 
representation for the ellipticity of the lensing cluster: 

£Ai689 = Y+J &2 ' a (12) 

and minimised £ in £i — £2 — $e parameter space. Fitting 
this profile shows that the data require an asphericity in 
the lensing mass distribution, which is also suggested both 
by the luminous mass distribution and by non-parametric 
mass reconstruction. 

For the NFW model in an Q = 1.0, A = 0.0 Universe, 
confidence contours (2A£ = 2 [£{ tt )—£ ( 7 r max )]) are shown 
in Fig- 0 - Confidence contours for the power-law model are 
shown in Fig.0. 

Fig-0 shows the enclosed projected surface mass Mj(< 
R) (in units of Mq) for each of the SIS, POW and NFW 
models, as a function of projected distance from the profile 
centre R (in kpc) plotted out to R = r 2 oo = 1.14 h~ l Mpc. 
Note the good agreement between the total masses esti¬ 
mated by fitting each of the models: the maximum differ¬ 
ence beween the models is ~ 20 % at the virial radius and 
smaller for smaller radii. 

3.2. Assessing the goodness of fit 

Just how good a fit is the best-fit NFW model, 7 r max ? 
Specifying g(tf) for 7 r max , the catalogue was “unlensed” by 
transforming the observed e,; to intrinsic e? under this lens 
model. Taking this reference catalogue, the phases of e® 
were randomised, while preserving 1 9i, and |ef |, to obtain a 
new catalogue. This catalogue was lensed and the value of 
({((imax), where the subscript R denotes randomisation, 
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Fig. 1. This figure shows the confidence contours for the 
NFW profile family. We show contours of constant 2A£ = 
2 [£( 7r)—£ (7r max )]. Starting with the innermost contour, 
the levels are plotted at 68.3, 90, 95.4 and 99-% confidence. 


was determined for this realisation. This randomisation 
of the reference catalogue, lensing and determination of 
fn{n max) was repeated 10000 times, to build up statistics 
of $ := 2 [£fl(7r max ) — £(7r max )]; the cumulative probability 
distribution P(< $) is shown in Fig^. The value £(7r max ) 
(fit to the real data set) lies well within the distribution 
of £_R(7r max ) for the random phase data sets, with 43.9 
(56.1)% of realisations being more (less) consistent with 
the model. 

We repeated the process for the power-law model, and 
in this case, 44.3 (55.7)% of randomised data sets are more 
(less) consistent with the model. 

3.3. The influence of the aperture size 

In order to test the influence of the size of the data field 
on the recovered NFW parameters, the value of 0 out was 
decreased from 15'.4 and 7r max obtained for these cata¬ 
logues. The largest deviation of 7r max for the smaller aper¬ 
tures compared with the full aperture is at most 10% until 
0 out ~ 8'.0. For instance, when 0 out = lO'.O the deviation 
of both r 2 oo and c from their values when 9 out = 15'.4 
is well under 1%. Below about 8'.0, the shear signal be¬ 
comes noisy, and changing the aperture slightly can have 
a significant effect on parameter recovery. Another effect 
of decreasing aperture size is to expand the confidence 
contours of the best-fit model; this has been discussed in 
detail in SKE. 

When 0 out = 15'.4 and 6% is increased from l'.O to 2'.0, 
the best-fit NFW model still lies within the original 68.3-% 
confidence contour. 



0 E (aremin) 

Fig. 2. Same as Fig. [I] for the power-law family. 

Table 2. This Table shows 7r max and the corresponding 
value of M 200 for the NFW profile, assuming different red- 
shifts for the background sources, z. An f i = 1.0, A = 0.0 
Universe is denoted by EdS, and an U = 0.3, A = 0.7 
Universe by A-flat. 


Cosmology 

z 

r* 2 oo 

( h -1 Mpc) 

c 

M 200 

(10 14 h,- 1 M 0 ) 

EdS 

0.5 

1.27 

5.2 

7.8 

EdS 

0.75 

1.18 

4.8 

6.3 

EdS 

1.0 

1.14 

4.7 

5.7 

EdS 

1.25 

1.12 

4.6 

5.4 

A-flat 

1.0 

1.29 

4.8 

5.95 


3.4. The impact of background source redshift and 
cosmology 

Concentrating on the NFW profile to examine the impact 
of our choice of background source redshift and cosmol¬ 
ogy, Table |] shows the values of 7r max obtained with vari¬ 
ous assumptions. To facilitate the comparison between the 
different scenarios, we also give values for M 2 qo, the mass 
contained within a sphere of radius r 2 oo and having mean 
density 200p c . Results are shown for an 9 = 1.0, A = 0.0 
cosmology, and for an U = 0.3, A = 0.7 cosmology with 
z = l. 

If the background sources are at z = 1.25 rather than 
at z = 1.0, M 2 oo would be about 5% lower, and if the 
sources were at 2 : = 0.75 the mass estimate would be about 
10% higher. Working in an (1 = 0.3, A = 0.7 cosmology, 
and keeping z = 1.0, gives an increase of about 5% in the 
best-fit value of M 2 oq. Again, we can compare this result 
with the value obtained by Clowe & Schneider assuming 
the same cosmology and z = 1.0: the values of r 2 oo and 
c in that work correspond to an M 2 oo value that is 2.5% 
higher than that obtained here. 












King, Clowe & Schneider: Parameterised models for the lensing cluster Abell 1689 


5 



1 Q13 I'll' . I I . .. I l .. I ... I l l l I .. I 

0 t 200 400 600 800 1000 \ 

9 in R (h- 1 Mpc) r 200 


Fig. 3. The behaviour of the enclosed surface mass Mj(< 
R ) as a function of projected distance R from the cluster 
centre, plotted out to r 2 oo = 1.14 /i -1 Mpc for the best- 
fit SIS (solid line), POW (dotted line) and NFW (dashed 
line) models. The inner fitting radius, 9 ln is indicated by a 
vertical dashed line (the outer fitting radius, 0 out , is larger 
than r 2 oo)- On the upper horizontal axis, the best-fit SIS 
Einstein radius from weak lensing is denoted by 9e 1 and 
the Einstein radius implied by the strong lensing arcs by 
9e (strong). 


4. Discussion and conclusions 

The mass profile of a cluster can be constrained by using 
the weakly lensed images of background galaxies; since 
lensing does not depend on the nature of the matter, 
it is an invaluable means to study the luminous and 
dark matter. One technique that harnasses the infoma- 
tion contained in the distorted ellipticities of the back¬ 
ground population is the shear likelihood method devel¬ 
oped in SKE and in KS. In this paper we have applied 
the method to WFI observations of the z,\ = 0.18 lensing 
cluster Abell 1689, to obtain best-fit SIS, SIE, NFW and 
power-law models. Up until now, we had rigorously tested 
our method on synthetic data sets, generated by lensing 
through analytic profiles and numerically simulated clus¬ 
ters. 

The maximum likelihood method developed in SKE 
and extended in KS differs from standard y 2 fitting [used 


Fig. 4. Goodness of fit results for the NFW model: the 
distribution of P(< $) for := 2 [f fi (7r max ) - £(7T max )] 
determined from 10000 randomisations of the background 
galaxy position angles. 

recently by Clowe & Schneider (2001) for example] in a 
few respects: 

— Each galaxy is treated individually rather than being 
binned. 

— The most likely model is found by finding the parame¬ 
ters that maximise the probabibility of the lensed ellip- 
ticity distribution. In y 2 fitting, the difference between 
the model and observed reduced shear is minimised, 
taking into account the observational errors. 

— The unlensed ellipticity distribution is assumed to be a 
gaussian in the maximum likelihood method, whereas 
the lensed ellipticity distribution is assumed to be a 
gaussian in y 2 fitting. 

Since a subset of the models considered here, namely 
the NFW and SIS models, were considered by Clowe & 
Schneider (2001), we can directly compare the results. 

The distorted images of background sources, at pro¬ 
jected distances between l'.O and 15'.4 (0.12/i” 1 Mpc < 
R < 1.8 /i _1 Mpc) from the centre of the cluster have been 
used to obtain the best-fit models. In an f l = 1.0, A = 0.0 
cosmology, the best-fit 1-paranreter SIS model has ctid = 
1028l|2kms -1 . This velocity dispersion is consistent with 
that of the X-ray study of Allen (1998), who investigated 
a cluster sample and obtained temperature profiles by ap¬ 
plying a deprojection technique to X-ray surface bright¬ 
ness profiles. Dynamical analyses of the cluster indicate 
that it is not a simple system: Teague, Carter & Gray 
(1990) found ctid = 2355kms _1 , and evidence for three 
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other groups of galaxies in its foreground. A more re¬ 
cent analysis by Girardi et al. (1997) finds that A1689 
is composed of three structures for which the resultant 
a id ~ 560km s -1 is much lower. Girardi et al. (1997) 
point out that their result, along with the foreground 
groups suggested by Teague et al. (1990), may indicate 
the presence of a filamentary structure along the line of 
sight. The strong lensing giant arcs imply an Einstein ra¬ 
dius 9e ~ 0'.75 (e.g. Miralda-Escude & Babul (1995)). 
Assuming the arcs are at z = 3, the corresponding SIS 
has ctid = 1360kms _1 . In fact, ctid is rather insensitive 
to z, decreasing by less than 2% even when z = 6. 

As discussed in Clowe & Schneider (2001), this data set 
gives lower limits for the mass of the lensing cluster, due 
to contamination of the galaxy catalogue with stars and 
foreground and cluster dwarf galaxies. Even with any re¬ 
alistic correction considered, they find that the weak lens¬ 
ing value can not be reconciled with that interpreted from 
strong lensing data. This points to the need for more spec¬ 
troscopic and photometric redshift information in both 
regimes, and a thorough combined analysis of the weak 
and strong lensing data. 

Formally, the best-fit 2-parameter model belongs to 
the NFW family (r 2 oo = l-14/i _1 Mpc, c = 4.7), but it 
cannot be distinguished from the other best-fit models at 
the 3 -ct level. In order to distinguish profiles at a satis¬ 
factory significance, so testing the CDM paradigm (that 
predicts NFW and similar cuspy profiles) stacking the sig¬ 
nals from a number of clusters is required. 

The best-fit general power-law has a slope q = 0.88, 
and 0e = 0'.3 (0.035 /i -1 Mpc). One might have expected a 
steeper slope, since most of the background galaxies sam¬ 
ple the profile at radii outside r s , where the slope is ex¬ 
pected to be steeper than isothermal. Our result indicates 
that galaxies close to 0; n are important in the determina¬ 
tion of the best-fit model. 

The confidence contours for the NFW result show that 
?’ 2 oo is better constrained than c, consistent with our study 
of NFW profiles in KS. Thus M 2 oo, which depends on r 2 oo, 
is also well constrained. When the power-law is fit, the 
confidence contours show a smaller fractional uncertainty 
in q than in 6e, which is in agreement with our finding in 
SKE. 

We also considered a 3-parameter SIE model; this 
model reveals that the projected (dark and luminous) sur¬ 
face mass density is not circularly symmetric. The ap¬ 
pearance of the cluster on optical images, and a direct 
non-parametric mass reconstruction also suggest deviation 
from circular symmetery. It is interesting that although 
the SIE is a 3-parameter model, its likelihood is in fact 
lower than that of the 2-parameter NFW model. This sug¬ 
gests that in order to fit the weak lensing data, any devi¬ 
ation from circular symmetry is not as important a driver 
as the deviation from a purely isothermal profile slope. 
Another factor coming into play may also be the “twist¬ 
ing” of the isodensity contours as a function of distance 
from the centre; this is evident when the fitting aperture 


is changed, and even directly on the non-parametric re¬ 
construction. 

By randomising the phases of the background sources, 
we tested that the NFW model and power-law models are 
indeed good fits to the real data. One point to note is that 
the galaxy catalogue must be unlensed before the phases 
are randomised and the catalogue is again lensed. If the 
catalogue is not unlensed, this results in a significant offset 
between the likelihood values from these randomisations 
compared with the maximum likelihood value from the 
original catalogue. Although only a small difference in the 
ellipticity dispersion is brought about by lensing, this is 
summed over 19400 galaxies when calculating l\ 

Note that changing the value assumed for the back¬ 
ground source redshift, z, only has a marginal effect on 
the best-fit NFW parameters: changing z by 0.25 results 
in an approximately 5-10% change in the M 200 value in¬ 
ferred from the best-fit values of r 2 00 and c. This indicates 
that our analysis is not too badly affected by the lack of 
redshifts for the background sources. However, we hope to 
obtain photometric redshifts for many of the background 
sources in order to obtain as accurate a result as possible. 
Assuming an Q = 0.3, A = 0.7 cosmology changes the 
mass estimate by « 5%. Again, this mass estimate is con¬ 
sistent (2.5% difference) with that of Clowe & Schneider 
( 2001 ). 

An important point to note is that the difference in 
the projected mass Ms of A1689 within R = r 2 00 for 
the best-fit SIS, power-law and NFW models is less than 
20%. This means that adopting a particular parameterised 
model does not lead to a severe error in the mass estimate. 

To conclude, we have applied the shear method dis¬ 
cussed in SKE and in KS to the spectacular lensing clus¬ 
ter Abell 1689. The algorithm requires only a few seconds 
of epu time to obtain best-fit parameterised models. Our 
results for the best-fit ctid for the SIS and Af 2 oo f° r the 
NFW profile are consistent with standard y 2 fitting; one 
main advantage of the maximum likelihood method is that 
galaxies need not be binned. We are very encouraged by 
the performance of the method, which will be applied to 
a larger sample of clusters in due course. 
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